scaling law in target-hunting processes 
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We study the hunting process for a target, in which the hunter tracks the goal by smelling odors 
it emits. The odor intensity is supposed to decrease with the distance it diffuses. The Monte Carlo 
experiment is carried out on a 2-dimensional square lattice. Having no idea of the location of the 
target, the hunter determines its moves only by random attempts in each direction. By sorting the 
searching time in each simulation and introducing a variable x to reflect the sequence of searching 
time, we obtain a curve with a wide plateau, indicating a most probable time of successfully finding 
out the target. The simulations reveal a scaling law for the searching time versus the distance to the 
position of the target. The scaling exponent depends on the sensitivity of the hunter. Our model 
may be a prototype in studying such the searching processes as various foods-foraging behavior of 
the wild animals. 
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In the past years, the diffusion-controlled reactions have been extensively studied through random-walk models. 
Such applications range from chemical processes, electronic scavenging and recombination, to electronic and vibra- 
tional energy transfer in condensed mediaQ, E&SSUHHIi. Man y 

works have been devoted to the target 
annihilation problem, in which randomly placed targets are annihilated by random walkers, and its dual of the trap- 
ping problem|lCl 1 1 1 1] . Other models treat hindered diffusion problems which involve random point obstacles [12| . In 
these models, the tracer moves from site to site on a lattice and falls into wells of various depth at the sites. The 
tracer does not know the depth of a well before it enters. Another possibility is a mountain model, in which all sites 
are at zero energy and the barriers are on the bonds joining the sites|l3t ll l| . Generally, random- walk models are 
ideally suited for computer simulations, a practical way to obtain results, since for the vast majority of cases no purely 
analytical method exists. 

In this work, we focus on another class of random-walk problems. We study the so-called target-hunting processes, 
which frequently occurs in biological systems, such as a shark searc hing for foods by smelling the blood in the ocean, 
or honeybees flying in the countryside to locate the foraging-nectars |l51 Il6l Il7| . or in metabolic processes such as cell 
motions and chemotaxis pll. fl9l l20l l2ll |22| . It can be viewed as the target-oriented problems, in which the hunters 
try to reach the targets by following some kind of behavior rules. In our model, an active hunter is trying to find out 
a target which emits a special kind of odor. The Monte Carlo simulations are carried out on a 2-dimensional square 
lattice. Since neither the distance nor the direction of the target is presumedly to be known, the searcher should 
determine its moves by random attempts in each direction, just like a snake turns its head from side to side to test 
the variation of the odor intensity. There are some chance for the hunter to move in the wrong direction because of 
randomness. Hence it is not a traditional biased random-walk. After sorting each searching process in time sequence, 
we obtain a curve with a wide plateau, indicating a most probable time of successfully finding out the target. By 
fitting the numerical results, we find a scaling law for the searching time on the distance to the position of the target. 
The scaling exponent is found to be dependent on the sensitivity of the hunter. We consider this scaling law scarcely 
happens in ordinary biased random walks. 

The game rules are as follows: The hunter at the origin O is trying to find out a target which emits a special kind 
of odor. Since the hunter have no way to know the location of the target, it randomly moves around its original 
position to test the variation of the odor intensity, zq is the present distance of the hunter to the target while z\ is 
the corresponding distance of the next attempted step. The Monte Carlo steps are implemented as: if {z( S Jz\) a > (, 
where the parameter a reflects the sensitivity of the hunter and £ is a random number, then the attempt is accepted. 
Otherwise it is refused. This rule implies that the intensity of the signal emitted by the target is inversely proportional 
to the distance of the hunter to the target. Other choices of the relation do not alter the result qualitatively. By this 
way, the hunter approaches the goal in a stochastic style. 

Fig.l displays a typical route of the hunter searching for the target on a regular lattice. When the hunter is far away 
from the goal, the ratio Zq/zi is close to 1. Most of the moving attempts are accepted, even the hunter walks in the 
wrong direction. The hunter appears to linger around for quite a while. Hence the motion of the hunter is nearly a 
Brownian random walk. As the goal is nearer, the ratio of zq/ z\ gradually approaches 0.5 and the probability of being 
refused for the hunter moving in the wrong direction increases. Hence the searching route seems more straightforward. 

Fig. 2 shows the searching time for each simulation for a distance of z = 31.4 and a = 6. It is understandable that 
the searching time are different for different stochastic processes. The distribution is not like a white noise. There are 
large fluctuations away from the most probable searching time. In Fig. 3, we plot the distribution of searching time 
for z = 31.4,65.6, 137.1, 188.4, respectively. It is seen that the distribution is not of the Poissonian form. The curve 
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has a very long time tail. Instead, a power relation is found for the maximal value of the distribution with the most 
probable time, 



V m ~i- d ™, (1) 



with d m — 1.05 for a = 6. 

There is another power relation between the distance of the target and the most probable searching-time, 

t P ~ z d ", (2) 

where d p = 1.77 for a — 6. 

In Fig. 4 we redistribute the data in Fig. 2 by sorting with increasing time. Fig. 4(a) is for various distances z from 
the origin, with z — 31.4,65.6,137.1,188.4 from bottom to top and a = 6. The horizonal axis is the sequence of 
searching time represented in percentage. A wide plateau is formed in the intermediate range. Fig. 4(b) shows that 
after proper displacement, all of the curves collapse into one, implying these curves be parallel to each other. Hence 
each curve can be described by a single function /(a;, a) plus a z-dependent function </>(z, a), 

lnt(x, z, a) = f[x, a) + (p(z, a) (3) 

with x the sequence of searching time represented in percentage. The function f(x) is of Arabic Ogive-like. It can be 
checked that 

lnt(^/zi ■ z 2 , x) = — (\nt(zi, x) + ln/j(z2, x). (4) 

Formula @ can be written as 

ln/j(x, z, a) = f{x, a) + 77(a) lnz. (5) 

In figure 5 we studied the dependence of the searching-time with respect to the sensitivity parameter a. Fig. 5(a) 
shows the curves for a = 4, 8, 12, 16, 20. After properly rescaling the curves in (a) by times Inrj with a coefficient or 
where f3 is determined below, all of the curves become parallel. From Fig. 5(b), we deduce or In rj = f{x) + (f>(z,a). 
By comparing with eq.JSJ), one gets 

lnt(a;, z, a) — oT^ f{x) + r/(a) lnz. (6) 

The index /3 can be derived by considering the dependence of the slope hi of the plateau in figure 5(a) on parameter 
a. There is a good linear relation between ln^2 and In a, as show in the inset of figure 6, 

lnfc 2 --/31na. (7) 

We measured (3 — 0.623. It is noteworthy that (3 is a constant independent of sensitivity parameter a. It results from 
the stochastic process. 

Finally, we try to found out the relation between 77 and a. Fig. 7(a) depicts the relation of hit versus lnz at x = 0.6 
for various parameter a. From figure 7(b), 

In 77(a) ~ —8a (8) 

with 5 = 0.01. 

Combining all the above factors, we consequently obtain a complete relation of the searching-time with respect to 
the distance as well as sensitivity parameter a, 

lnt(x,z,a) = a- p f{x)+c Q e- 6a \nz, (9) 



t[x,z,a) =e f{x)/a " -z d . (10) 
We find that there is a generalized scaling-law between t and z with exponent 

d = c Q e- Sa . (11) 
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From formula J3J, a = 6, d = d p = 1.77. We get Co = 1.88. It shows that the power-law exponent is a-dependent. As 
a increases from zero to infinity, the exponent decreases from 1.88 to zero. In Ea. (|10[) . the contributions of variable 
x, which sorts the searching time in each simulation, are completely merged into a prefactor and the scaling exponent 
is x-independent. It should be noted that it is a functional relation between the searching time and the distance in 
such stochastic processes. 

In summery, we introduced a variable x to denote the sequence of searching time. We plot a curve with a wide 
plateau, indicating a most probable time of successfully finding out the goal. In stead of calculating the mean square 
root, we introduce a sort parameter x to figure out an analytical expression. The simulations reveal a scaling law 
for the searching time versus the distance to the position of the target. The scaling exponent is dependent on the 
sensitivity of the hunter. We believe that our treatment of the statistical data may be useful in other cases. The 
existence of the scaling law may have implications with the possibility for the hunter to walk in a wrong direction 
or stay at the same place for quite a while. It scarcely happens in an ordinary biased random walk. We point out 
that the results are valid not only on the square lattice, but also for continuous moving (with fixed step length) in 
the two-dimensional plane. However, the explicit form of function f(x) is still lack. It is also desirable to deduce an 
analytical express of Eq. I|10(l from the first principle of statistics. 

We suggest that the scaling law in the hunting process may be an additional behavior rule in the foods-foraging 
processes of wild animals, which has not caught much attention. In turn, verifications of the law from direct obser- 
vations by zoologists or entomologists are also expected. Our target-oriented model may be a prototype in studying 
the foods-foraging processes in wildlife as well as in other searching games. 
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Figure Captions 

Figure 1 A typical hunting route on a square lattice. The start point is at the origin and the target is at (25, 19). 
Figure 2 Searching time for 10000 simulations. The original distance to the target is z = 31.4 and the sensitivity 
parameter a — 6. 

Figure 3 Searching time distribution for distances z = 31.4,65.6, 137.1, 188.4, respectively, a = 6. 

Figure 4 Time-sorted curves for 10000 simulations. The horizontal axis x is the sequence of searching time repre- 
sented in percentage and the vertical axis is logarithmic time, a — 6. (a) is for various distances z from the origin, 
with z = 31.4,65.6, 137.1, 188.4 from bottom to top. Evidently, the curves consist of three parts, and a wide linear 
region is formed. These curves is parallel to each other, (b) shows that after proper displacement, all curves collapse 
into one. 

Figure 5 Time-sorted curves for 10000 simulations for various parameter values a. The horizontal axis x is the 
sequence of searching time represented in percentage and the vertical axis is logarithmic time. The original distance 
is fixed at z = 31.4. (a) is for a = 4,8,12,16,20 from top to bottom, (b) shows the rescaled curves of (a) for 
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a = 4,8, 12, 16, 20 from bottom to top. These curves are parallel to each other. 

Figure 6 The dependence of the slope of the plateau in figure 4 on parameter a. The distance z = 31.4. Inset: a 
linear relation of the logarithmic slope ki with logarithmic a. 

Figure 7 (a) A plot of hit versus lnz at x = 0.6 for various parameter a. (b) A linear relation of In 77(a) with 
respect to a. 
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